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Abstract 

In studies of interfaces with dynamic chemical composition, bulk and interfacial quantities are often 
coupled via surface conservation laws of excess surface quantities. While this approach is easily justified 
for microscopically sharp interfaces, its applicability in the context of microscopically diffuse interfaces 
is less theoretically well-established. Furthermore, surface conservation laws (and interfacial models 
in general) are often derived phenomenologically rather than systematically. In this article, we first 
provide a mathematically rigorous justification for surface conservation laws at diffuse interfaces based 
on an asymptotic analysis of transport processes in the boundary layer and derive general formulae for 
the surface and normal fluxes that appear in surface conservation laws. Next, we use non-equilibrium 
thermodynamics to formulate surface conservation laws in terms of chemical potentials and provide 
a method for systematically deriving the structure of the interfacial layer. Finally, we derive surface 
conservation laws for a few examples from diffusive and electrochemical transport. 

Keywords: surface conservation laws; transport theory; interface structure; asymptotic analysis 

1 Introduction 

It is well-known that the chemical composition of an interface plays an important role in its dynamics 
(both its motion and the evolution of its chemical composition) . A few interesting and important examples 
include the effect of surfactants on surface tension [TJ [H [3J , the charging dynamics of electrochemical double 
layers H [5J H [D [SJ , electrokinetic phenomena PlIHlinilllllllllllSlIISlIIlIISlIIllini and soap film 
dynamics [211 122j . Unfortunately, theoretical modeling of the interaction between chemical composition 
and interfacial dynamics can be quite challenging because it involves coupling dynamics at the macroscopic 
scale away from the interface with a microscopic model for the interfacial region. Sometimes it may not 
even be clear which microscopic model is appropriate for a particular problem. Many successful approaches 
based on the notion of an excess surface concentration HI] have been developed and applied to a wide- 
range of problems. While this approach is well-founded for microscopically sharp interfaces, the theoretical 
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foundations for this methodology in the context of microscopically diffuse interfaces (especially for dynamic 
interfaces) do not appear to have been thoroughly explored. 

For microscopically diffuse interfaces (i.e., interfaces where the chemical composition smoothly varies 
across the interface) , the notion of a surface excess concentration and dynamics of surface excess quantities 
can be rigorously justified via asymptotic analysis. In this article, we present a general procedure for deriving 
surface conservation laws (including derivation of the structure of the interfacial layer) that couple dynamics 
in the bulk with dynamics of the interfacial region and demonstrate the application of our formulation to 
example problems in diffusive and electrochemical transport. Our results are only valid in "sharp interface" 
limit (i.e., when the distance over which the volume density of the diffuse species varies is small relative to 
macroscopic length scales of the problem). Fortunately, this limit is the only one that is physically meaningful 
- if an interface is too diffuse, it may not be appropriate to treat it as an interface in the first place. 

Previous work involving surface conservation laws focuses primarily on deriving the appropriate surface 
conservation law in the context of specific problems [TQl HI [7] . In contrast, our derivation is much 
more general and applicable to a wide range of transport problems. We only require that, near the surface, 
the flux for the transport process scales in the same manner as a linear combination of gradients of field 
variables. Fortunately, these types of fluxes are very common in transport theory [25] . We emphasize that the 
formulation of surface conservation laws we present is not based on ad-hoc physical arguments; rather, they 
are a direct consequence of the asymptotic analysis. Moreover, the generality of our formulation allows us to 
derive surface conservation laws for nonlinear transport processes when the bulk is far from equilibrium (e.g., 
electrochemical transport of around metallic particles at large applied fields) and for systems with complex 
interfacial structure. To the authors' knowledge, surface conservation laws for these types of problems have 
not been previously derived. 

1.1 Dimensionless Formulation of Equations 

To facilitate the mathematical analysis presented in this article, it is convenient to examine all equations in 
nondimensional form. To avoid confusion when reintroducing dimensions, let us take a moment to fix the 
physical scales used in the nondimensionalization process. At the macroscopic scale, the spatial coordinates, 
(X, Y, Z), are scaled by the the characteristic size of the system, L. At the microscopic scale, the tangential, 
(x,y), and normal, z, spatial coordinates are scaled by L and S, respectively, where S is the characteristic 
thickness of the interfacial layer. Time, t, is scaled by the bulk diffusion time scale, r = L 2 /D. Concen- 
trations, Cj, are nondimensionalizcd by taking the reference concentration to be the maximum physically 
realizable concentration, c max . Bulk, Fj, and surface fluxes, J^, are scaled using c max L/r and c max L 2 /r, 
respectively. The electric potential, <fi, is scaled by the thermal voltage fc^T/e, where fee is Boltzmann's 
constant and e is the absolute charge of an electron. Finally, all thermodynamic energy variables (e.g., U 
and F), and entropy, S, are scaled using fcgT (the thermal energy) and fcs, respectively. The dimensional 
form of any of the equations in the remainder of this article can be obtained by first replacing each dimen- 
sionless variable with its associated dimensional variable divided by the appropriate characteristic scale and 
then multiplying the entire equation by an appropriate physical scale to restore units to the equation 1 . For 
instance, if the dimensionless concentration, c, appears in a dimensionless equation, it should be replaced by 
C I c m ax when re-dimensionalizing the equation. 



2 



Table 1: Summary of physical scales used to nondimensionalize equations. 
Physical Quantity Variable(s) Physical Scale 



spatial coordinate 


X, Y, Z, x, y 


L (macroscopic length scale) 


spatial coordinate 


z 


6 (microscopic length scale) 
t = L 2 /D (bulk diffusion time scale) 


time 


t. 


concentration 


Ci 


Cmax (maximum physically realizable concentration) 2 


bulk flux 


F- P 




surface flux 






electric potential 


<a 


ksT/e (thermal voltage) 


energy 


U,F 


ksT (thermal energy) 


entropy 


S 


k B 



2 Excess Surface Concentrations 

When studying chemically active interfaces, the state of the interface is commonly described at macroscopic 
scales by specifying the excess surface concentrations of all chemical species. Intuitively, excess surface 
concentrations are defined as the amount of material per unit area of surface after the material in the "bulk" 
has been removed from the interfacial region. However, this definition is physically and mathematically 
ambiguous - in the region near the interface, how is one to distinguish between material that is part of the 
bulk and material that is part of the interface? The ambiguity of this definition was recognized long ago by 
Gibbs [5Sl[571[2Hj and is typically dealt with by arbitrarily selecting a concentration, C* , at some point within 
the interface (when viewed at a microscopic length scale) to be the reference "bulk" concentration .23, 2~Hl2T)]. 
Any deviation of the concentration near the interface from C* is treated as a contribution to the excess surface 
concentration. In terms of this reference concentration, the excess surface concentration is defined as the 
integral of (C — C*) over the thickness of the interface. 

While somewhat inelegant from a theoretical perspective, this formulation of excess surface concentration 
has been used quite successfully for studying interfaces of equilibrium systems. For systems in thermal 
equilibrium, the excess surface concentrations can be related to other thermodynamic variables, such as 
bulk concentrations and other excess surface concentrations via adsorption isotherms [30l [IJ [23l EH [31] . 
However, when the dynamics of surface species is important {e.g., fast adsorption-desorption kinetics, non- 
negligible surface-diffusion) , isotherm models need to be replaced by surface conservation laws for the excess 
surface concentrations which provide dynamic coupling between bulk concentrations and excess surface 
concentrations : 

^- = -V s -J s + Jn- (1) 

ot 

Here T, J s and J n are the surface excess concentration, the surface flux and the normal flux, V s denotes a 
surface derivative, and the sign on the normal flux is chosen to be positive when the flux is into the boundary 
layer (see Figurc[T]). In this situation, the thermodynamic definition of excess surface concentrations becomes 
unsatisfactory because all of the variables in ([T]) depend on implicitly on the choice of C* whose governing 
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equation and relationship to other macroscopic variables may be difficult to derive. 




Figure 1: Schematic diagram of fluxes at an interface. Note that the surface excess 
concentration T and the surface fluxes J s are defined only on the interface itself. 

Fortunately, for many problems, there is a natural way to define C* in terms of bulk concentrations 
(which almost always have relatively straightforward governing equations). Define C* to be the limit of the 
corresponding bulk concentration, C(x), as x approaches the interface. Unlike thermodynamic definitions 
of the excess concentration which arc defined in terms of a reference concentration chosen at microscopic 
length scales and suffer from the difficulty of relating that reference concentration to macroscopic variables, 
the choice of C* just described is directly related to macroscopic variables because it is explicitly defined in 
terms of those variables. Note that C* is not necessarily a constant; it may still be a function of position 
along the surface of the interface. 

3 Surface Conservation Laws 

Equation (fTJ) seems physically intuitive: the time rate of change in the surface concentration results from 
a combination of surface diffusion and flux from the bulk. However, it is important to remember that the 
equation describes the evolution of excess concentrations, not absolute concentrations. Thus, 3 S and J n must 
be carefully defined so that they contribute solely to changes in the excess concentration, not to changes in 
reference concentration C*. 

For microscopically sharp interfaces, such as the monolayer interfaces that arise in problems involving 
surfactants at liquid-gas interfaces [H [2j [3] or the compact layer in electrochemical systems [30l [24] , J s and 
J n are very simply defined because the bulk truly extends all the way to the interface. As a result, the excess 
surface concentration is the amount of material per unit surface area that resides precisely on the interface. 
In other words, the excess surface concentration is the absolute surface concentration. Thus, the surface 
can truly be thought of as a distinct phase, and a conservation law argument based on a balance of fluxes 
into a "control volume" on the interface is mathematically valid. As a result, J n is the flux of bulk material 
normal to the surface and J s is the flux within the surface itself driven by transport processes intrinsic to 
the interface (e.g., surface diffusion in an lipid bilayer). 

That surface conservation laws hold for microscopically sharp interfaces is not surprising. What is 
interesting is that J s and J„ can be defined in such a way that surface conservation laws also hold for 
microscopically diffuse interfaces (see Figure [2]). For instance, in their studies of the evolution of surface 
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Figure 2: Schematic diagram of a microscopically diffuse interface. The dashed line 
represents the outer "edge" of the boundary layer, which is not at a mathematically 
well-defined location. Note that in addition to the fluxes within the boundary 
layer, it is possible for the interface to possess fluxes that physically reside on the 
microscopic interface itself. For some systems, these "microscopic" surface fluxes 
may provide a significant contribution to the total surface flux in the boundary 
layer. 



excess concentrations of electrolyte around colloid particles [lOl IH H] , Dukhin, Deryagin, and Shilov made 
use of surface conservation laws for excess ion concentration to couple ion transport in the bulk to the 
behavior of the electrical double layer. A key observation about microscopically diffuse interfaces is that 
surface transport within the boundary layer near the interface is driven by the same transport processes that 
occur in the bulk. In contrast, at microscopically sharp interfaces, surface transport could potentially be 
completely different from bulk transport because the underlying physical processes driving surface transport 
(which are related to the specific nature of the interface) might be fundamentally different. 

Unfortunately, the derivation of surface conservation laws for microscopically sharp interfaces is not valid 
for diffuse interfaces because the notion of a surface "control volume" is no longer well-defined. The main 
problem is the lack of a distinct separation between the bulk and the interfacial region. Rather, there is a 
thin region near the interface over which concentrations vary rapidly. Intuitively, what we would like to do 
is define surface quantities by directly integrating over the boundary layer. This procedure, however, must 
be carried out carefully in order to obtain physically meaningful results. The basic idea in deriving surface 
conservation laws, which we shall elaborate on in the next section, is that all integrations over the boundary 
layer should involve only excess quantities; integration of absolute quantities leads to divergent results. 

Before deriving the general form for 3 S and J n involved in surface conservation laws for microscopically 
diffuse interfaces, we mention that {T]) neglects contributions due to motion and deformation of the interface 
(e.g., terms involving convection in the normal direction and surface dilation [II, 32, 33 ). Since our analysis 
sheds no new insight for these terms, we shall ignore them to simplify the discussion. 

3.1 Derivation of Equation (CD) for Microscopically Diffuse Interfaces 

To derive the form of the fluxes for surface conservation laws at microscopically diffuse interfaces, our goal 
is to identify the appropriate definitions for J s and J n (in terms of macroscopic variables) that will allow us 
to write surface conservation laws of the form ([1]). We begin by writing the conservation laws that govern 
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the time evolution of a chemical species in the bulk: 



dc 
di 



-V • F 



(2) 



where c is the concentration and F is the flux. For our derivation, we assume that there is no flux of 
material through the interface 3 itself so that dF/dn = and that the spatial coordinates in ((2]) have 
been nondimensionalized using a characteristic macroscopic length scale. In these units, the thickness of 
the boundary layer over which c varies rapidly is 0(e), where e = 8/L is the ratio of the interfacial and 
macroscopic length scales. 

The "sharp interface" limit is defined as the limit e — > and is best studied via asymptotic analysis [34| [35] 
using |j2]) for the outer equations and deriving the inner equations by appropriately rescaling the spatial 
coordinates. In this discussion, we shall distinguish inner and outer variables by using the hat (") and 
bar accents (~), respectively. Also, we shall use (X,Y,Z) (upper case variables) and (x,y,z) (lower case 
variables) to represent the outer and inner spatial coordinates, respectively. Since the focus will be on the 
region immediately neighboring the surface, we shall take (x, y) to be Cartesian coordinates tangential to 
the interface and z to the coordinate normal to the interface. 

Rescaling ([2]) using the inner coordinates (x, y, ez) = (X,Y,Z), we find that the governing equations in 
the boundary layer are 



where the subscripts s and n indicate tangential and normal components of the flux and divergence operator, 
respectively. Note that in ([3]), we have implicitly assumed that in changing to the inner coordinates, the 
normal derivative of the normal flux F n picks up a factor of 1/e because 



This relationship between the inner and outer normal fluxes is a direct consequence of our requirement that, 
near the interface, F scales in the same manner as a linear combination of gradients of macroscopic field 
variables: 



where fi are arbitrary functions of c\ , C2, . . . , c n . For fluxes of the form ([5]) , the factor of e in iQJ arises 
from the fact that the flux itself involves derivatives of space. We emphasize, however, that for our general 
analysis, we only require that fluxes satisfy (j4|) without regard for the origin of this relationship. 

In deriving ([1]), the intuitive idea of just integrating ([3]) from z — to z — oo is inadequate because it 
leads to divergent integrals which are physically and mathematically meaningless. The key idea to keep in 
mind when asymptotically integrating over boundary layers is that only excess concentrations are integrable, 
absolute concentrations are not. 

Using this basic principle, we can systematically derive (pj by starting with the excess (volume) concen- 
tration in the boundary layer: 



F n {x,y,z) = eF n {X,Y,Z). 



(4) 




(5) 



7(x,y,z,i) 



c(x,y,z,t) - c*(x,y,t), 



(6) 
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where c*(x,y,t) = lim^-^o c(X, Y, Z, t). Notice that, as discussed in Section ® we have chosen the reference 
concentration as limit of the bulk concentration as the interface is approached. Taking the time derivative 
of this equation and using inner and outer evolution equations, ^ and ([3]), we find that 

^)-(-v.-:-S§). 

where F* = lim^^o F. To obtain an equation for the surface excess concentration T, we would like to 
integrate this equation over the entire boundary layer. Unfortunately, it is not possible to just integrate over 
the entire range of the inner variable because the integral of the last term on the right hand side is divergent: 

lim dz= ( Um r dz. (8) 




z^o 8Z J \z^o dZ r 

However, because the boundary layer has an 0(e) width and all of the integrands are 0(1), we expect that 
all of the integrals should be 0(e) quantities. 

The problem with the intuitive approach is that it makes the mistake of equating the asymptotic limit 
e — > with the spatial limit z — > oo. Realizing this subtle distinction (which is safe to neglect for many 
asymptotic analyses), we can reformulate the integration over the boundary layer as the limit of a sequence 
of integrals over finite intervals, which tends to the entire half-space [0, oo) as e — > 0. In choosing the domain 
of integration, we want to be sure to capture the entire boundary layer so that the notion of the total 
surface excess concentration is physically meaningful. In addition, we want the region of integration at the 
macroscopic length scale to go to zero as e — > so that we are truly integrating over only the boundary layer. 
We can simultaneously achieve both of these goals by taking the region of integration at the macroscopic 
length scale (i.e., where Z is the variable in the normal direction) to be [0, a] where a = 0(e p ) with < p < 1. 
Since the width of the boundary layer is O(e), this choice of a ensures that the integration region completely 
covers the boundary layer but tends to in the asymptotic limit. 

Even with this choice of integration region, we still must exercise care to make sure that all integrands are 
of excess quantities so that integrations over the inner coordinate are convergent. This restriction suggests 
that we rewrite (0 as 

d-f /- \ 1 dF n 8F„ 

— = -V s • I F s — F* ) — — — h km (9) 

dt V / e 2 dz z^o dZ 

Integrating this equation over the boundary layer (at the macroscopic length scale), we obtain 
Changing from the outer to the inner coordinate for all of the integrals except the last term yields 



e 



7 dz = -e V,- (F.,-F!W-- / 4r^dz + lim ^? ) dZ 



dt 

Jo v ' s / e J dz J \z^o dZ 



-e 



a/ * , 1* , OF 



(i 



V.. (F s -F*) dz- -F n (x,y, a/e) +a I Um ) (11) 



Note that to move from the first to the second line in the above equations, we have explicitly integrated the 
normal derivative and applied tke no flux boundary condition, = 0) = 0. Expanding tke last term in 
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(fTTj) using a Taylor series, we find that 







^-dz = -e / V s - (F s -F*)<te- -F n (x,y,a/e) 



c 



+ F n (X, Y, a) - F n (X, Y, 0) + 0(a 2 ). (12) 
Finally, recalling that F n (x, y, a/e) = eF n (X, Y, a), the above equation can be simplified to 

-j+dz = -ej^ V,- (F s -F:jrfz-F„(X,r,0) + O(a 2 ). (13) 

By choosing 1/2 < p in the definition of a, we find that the 0(a 2 ) term becomes negligible compared to 
the remaining terms in the e — > limit so that the leading order asymptotic equation describing surface 
concentration evolution satisfies: 

^dz - -ej V.- (F s -F;)dz-F n (X,y,0) (14) 

Thus, by substituting the definition for the surface excess concentration 



we can make the identifications 



and 



jdz, (15) 



dz (16) 



J n = F(X,Y,0)-n = -F* n (17) 

to arrive at the surface conservation law ([T]). The sign difference between J n and F n is merely a byproduct 
of the choice of orientation for the local coordinate system in our analysis. As mentioned earlier, the sign 
convention for the normal flux is that J n be positive when the direction of the flux is into the boundary 
layer. 

It is worth mentioning that the presence of the e in the time dependent term and the surface flux term 
in (I14p indicates that the relative importance of these terms relative to bulk transport (i.e., the normal flux 
term) may depend on the choice of time-scales and the magnitude of surface transport. This observation is 
elaborated upon in the examples discussed below. 



4 Formulation in Terms of Chemical Potentials 

In general, it best to express driving forces for fluxes in terms of gradients of chemical potentials [313 [36] : 

Fj = -Ci b^LijVfij] +ua, (18) 

where the Lij are mobility coefficients 4 that relate the drift velocity of species i to the gradient of the chemical 
potential of species j and u is the background velocity that contributes to advective transport. Note that in 
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order to satisfy 0$, we require that the normal component of the background velocity, u n , at the interface 
must vanish. 

Substituting this expression into P^|) and (JTTJ) and rearranging a bit, we obtain 



(19) 



(20) 



where we have imposed u n = 0, u s is the tangential component of the background velocity, and the superscript 
i indicates on the 3 S and J n indicates that these fluxes are for the i-th species. 

The expression for J* can be simplified because the last integral term is equal to zero. As Q shows, 
a rescaling of the conservation law associated with (|T8|) to the inner coordinate yields the leading order 
equation: 



= 0. 



(21) 



Integrating this equation and using asymptotic matching for the boundary conditions as z — > oo, we find 
that at leading order, 

^Z'-'il'j ^I'-^l'j ( 22 ) 

3 3 

in the inner layer. That is, special linear combinations of the chemical potentials involving the mobility 
coefficients are constant in the normal direction and slowly varying in the tangential direction. Equation 
(I2ip has the following physical interpretation: in the sharp interface limit {i.e., e — > 0), the interfacial layer 
is always in quasi-equilibrium at the bulk diffusion time scale, r = L 2 /D. It is important to note that at 
faster time scales, which may be appropriate for systems with external forcing, this result may break down 
because the time derivative term could balance the normal flux term when ^ is rescaled to the faster time 
scale. 

Using this observation, J* becomes 



IX 



£i(u s — u*)dz. 



The surface conservation law for transport follows directly from these results: 



9£i 
dt 



= V,- 



Ci{u 8 — u*)dz 



(23) 



(24) 
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It is worth mentioning that the surface transport term (the term involving the surface divergence) does not 
always contribute to the leading order surface conservation law. Whether the surface transport term must be 
retained at leading order depends on the magnitudes of I\ (which depends implicitly on bulk field variables) 
and the tangential component of bulk chemical potential gradients. Interestingly, when the surface transport 
term is significant, the surface conservation law depends explicitly on the small parameter e (through 

4.1 Theoretical Derivation of Interfacial Structure 

While surface conservation laws may be derived by substituting arbitrary models for the structure of the 
boundary layer into (JTSJ) , (|16[) , (JTTJ) , the quasi-equilibrium nature of the interfacial layer allows us to the- 
oretically derive the boundary layer structure from the chemical potential using (|22|) . Starting with a free 
energy for the system that incorporates the physical effects we wish to include, chemical potentials may be 
easily computed as the functional derivatives of the free energy with respect to the concentrations of the 
individual species 5 : 

Hi = SF/Sc tl (25) 

where F is the total free energy of the system. Thus, formulating the transport equations in terms of chemical 
potentials yields a systematic method for deriving boundary layer structure and surface conservation laws 
directly from fundamental physical principles without resorting to ad-hoc models for the interfacial layer. 
We shall demonstrate this general approach for several example problems in the next few sections. 

5 Applications to Neutral Solutes 

In this section, we present surface conservation laws for neutral solutions. Because we have formulated surface 
conservation laws in terms of chemical potentials, derivation of the integrated fluxes for these systems is a 
straightforward application of the formulae in Section 01 For simplicity, we shall assume that fluid flow is 
absent. 



5.1 Dilute Solutions 

For dilute solutions of a single neutral species, we can write the free energy, F, for the system as 

F = U-TS= f c$ + c(logc- 1) dx, (26) 

where U, T, and S are the internal energy, the absolute temperature (scaled to 1 in dimensionless units), 
and entropy of the system, respectively, c is the concentration of solute, and $ is the energy of interaction 
between solute particles and the surface (per unit concentration of solute particles). Note that the free energy 
density per unit concentration of solute in (|26p is determined up to an arbitrary additive constant, which 
we take to be —1 {—kT in dimensional form) to simplify the expression for the chemical potential (derived 
below). The interaction energy, $(z), between solute particles and the surface extends over a distance e 
(typically a few molecular diameters) and can account for hydrophobic interactions, polarization effects, etc. 
$(z) essentially gives the affinity of the solute particles for the surface. Note that in the dilute solution 
approximation, solute particles do not directly interact with each other or with solvent particles. 
Computing the functional derivative of (|26p . we obtain an expression for the chemical potential 

/i = logc + $. (27) 
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As we showed in Section [H in the interfacial layer, p, is equal to fi* in the normal direction, which leads to 
a Boltzmann distribution for the concentration: 



c(z) = exp (fj,* - $(z)) . (28) 

We can rewrite this expression in terms of the bulk concentration just outside of the interfacial layer by 
observing that c* = e M from asymptotic matching and the fact that the solute surface interaction decays 
far away from the surface. Thus, we find that 

6(«) = c*e-*W. (29) 

Note that this is exactly the boundary layer profile obtained by Anderson et al. using the assumption that 
the interfacial layer is in equilibrium [37 . 

The surface excess concentration of solute, T, is the integral of the excess concentration over the boundary 
layer: 

r = e y (c(z) - c) dz = e J c* (e-* (z) - lj dz. (30) 

After the integral in this expression has been explicitly evaluated, the surface conservation law for dilute 
solutions of neutral solutes is derived by substituting the surface excess concentration (|3"U1) and the chemical 
potential (|2T|) into which leads to the following transport equation for the surface excess concentration: 

(31) 

where L is the nondimensional mobility coefficient for solute particles. It is interesting to observe that the 
driving force for surface transport and the source term for transport between the bulk and interfacial regions 
depend only on the concentration in the bulk. This result is generally true for transport problems where 
fluid flow is negligible. 



dt.- L 



v s -(r v s io gc *)-^ 

an 



5.2 Concentrated Solutions 

For concentrated solutions of a single neutral species, the size of the individual solute particles cannot be 
neglected and steric effects must be accounted for. Concentrated solution theory provides a model which 
qualitatively accounts for these effects. Bulk solutions are considered concentrated when the dimensionless 
parameter v — a 3 c max , which represents the bulk volume fraction of solute, is not small. In this case, we 
must use concentrated solution theory to describe bulk transport. However, even when v 1, concentrated 
solution theory may be necessary if the surface interaction energy, $, is strong enough to condense particles 
near the surface. As for dilute solutions, we begin by writing the free energy for the system using ideal 
solution theory |38j : 

F = U-TS = J c$ + [c log c+(l-c) log (1 - c)) dx, (32) 

where the last term in F accounts for then entropy of the solvent (i.e., steric effects) and we have assumed 
that solute and solvent molecules have the same size. Taking the functional derivative of this expression, we 
find that the chemical potential for concentrated solutions is given by 



ii = log ( — ) + $. (33) 
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Note that unlike dilute solutions, there is a maximum concentration for the solute that arises from the sharp 
increase in the chemical potential at high solute concentrations (low solvent concentrations) caused by the 
solvent's entropic contribution to the free energy. 

As before, we can derive the concentration profile of solute atoms by using the fact that, within the 
interfacial layer, the chemical potential is constant in the direction normal to the surface: 

c{z) = exp 0* - $(») [1 + cxp (jj,* - $(z))] _1 . (34) 

For concentrated solutions, the simple Boltzmann distribution for the quasi-equilibrium concentration profile 
in the normal direction is replaced by a Fermi-Dirac-like distribution with a maximum concentration of 1 
{cmax = 1/a 3 in dimensional form). Note that expression of c(z) in terms of the bulk concentration, c*, 
outside of the interfacial layer does not lead to a significant simplification of (j3"4")) . The surface excess 
concentration is found by integrating (|34|) over the interfacial layer. 

To derive the surface conservation law for concentrated neutral solutions, we substitute (|33[) and the 
expression for the surface excess concentration into the general formula for the surface conservation law in 
terms of chemical potentials (f2~4"|) . 



6 Application to Electrolytes 

In electrochemical systems, the solute particles are electrically charged, so transport processes are affected 
by electric fields (both self-generated and externally applied). When modeling electrochemical systems, it is 
common to model transport in the bulk separately from the response of the double layer at surfaces. Surface 
conservation laws justify and generalize the usual phenomenological approaches for coupling bulk and double 
layer dynamics. In this section, we discuss surface conservation laws for dilute and concentrated electrolyte 
solutions. As we shall see, the key difference between these two cases is the microscopic model that must be 
used to describe the electrical double layer. 

For electrolytes, transport is driven by gradients in the electrochemical potential, so the nondimension- 
alized conservation laws for ionic species are given by 

dc ' 

— = V • (c,V/i t + c A v) , (35) 

where Cj are the concentration of ion i, ^ is the electrochemical potential for species i, v is the fluid velocity, 
and we have assumed that all of the ions have the same diffusivity [31] . Note that the flux in this equation 
is just a special case of ([T8|) where Ly is an identity matrix 6 . While fluid flow plays an important role in 
many systems of current interest {e.g., electrokinetic microfluidic pumps and mixers |13 [ I14 | [T9] , 2(1). we 
shall neglect it to keep the following discussion simple. In this situation, the surface conservation law 
becomes 

§ = v,. ( r ( v,„.)-c.M. m 

Before considering ion transport in various special cases, we mention that for electrochemical systems, 
the small parameter e is the ratio of the Debye screening length, Ad, to the characteristic system size. In 
this article, the Debye length will be defined by 



" - 1/5*5- (37) 
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where e s is the dielectric constant for the electrolyte and C is the average concentration of neutral salt. 
Alternative choices for the Debye length typically differ from this definition by an O(l) multiplicative con- 
stant. 



6.1 Dilute Electrolytes 

The free energy for general dilute electrolyte solutions is similar to the free energy for dilute neutral solutions: 



F = U-TS= / ^QOogQ-lH^Tz^-e 2 ^! 2 )dx. (38) 



In this expression, the first term in the integrand is the entropic contribution from the solute particles, the 
second term accounts for the interaction energy between the charged particles and the electric potential, and 
the last term is the energy of the electric field. Note that e is defined to be consistent with the definition in 
([37]). 

Taking the functional derivative of the free energy with respect to Cj, we find that the electrochemical 
potential of the i-th ionic species takes the simple form 

[ii = logQ + Zi4>. (39) 

Using (|39[) for the electrochemical potential, (|3"5]) reduces to the commonly used Nernst-Planck equations [3D1 
[39]: 

= V ■ (Vc, + z lCl V4>) (40) 

Within the double layer, the electrochemical potential is constant in the normal direction, so the concentra- 
tion profile possesses a Boltzmann distribution: 

h{z) = c*e-^ z \ (41) 

where ip( z ) = 0( z ) — 

For the case of a symmetric, binary electrolyte, (|4Tj) leads to the commonly used Gouy-Chapman-Stern 
(GCS) model for the double layer (30j [23j [31] [40] . It is well known that the GCS model can be theoretically 
justified by noting that the double layer is always in quasi-equilibrium in the sharp interface limit [41] . 
However, the quasi-equilibrium structure of the double layer is not at all surprising in light of our discussion 
in Section[?]- it follows directly from the fact that the electrical double layer is a sharp interfacial layer. 

For symmetric, binary electrolytes, we can derive explicit expressions for the excess surface concentrations. 
In this case, the excess concentration of each of the two ionic species is given by 



7± = c ± 



(e Tz +$ - l) (42) 



where c* = (c5_ + cl)/2 is the average concentration in the bulk and z + is the charge number for the positive 
ion. Following [42] . this quantity is straightforward to integrate by changing the variable of integration in 
([T5]) to V> and using the fact that 
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for the GCS model. Carrying out the integration, we find that the excess surface concentration of species i 
is given by 

r± = ^H. (W/ 2 - i) (44) 



where C, — (f>(0) — 4>* is the zeta-potential [UJ 1301 131] across the diffuse part of the double layer. Using this 
expression for r± in (|36p . the surface conservation law for symmetric, binary electrolytes is fully specified 
in terms of bulk field variables and boundary conditions, and the structure of the boundary layer has been 
completely integrated out. 

Because the charge density and neutral salt concentration are both important components of the response 
of an electrochemical system, it is interesting to derive surface conservation laws for these quantities [H 02] . 
Toward this end, we define eq and ew to be the surface charge density and surface excess neutral salt 
concentration, respectively 7 : 



lf( T+ - 7 ,^^ + -r,=-^-(^) <«> 

™ = !^(7++7-)<fe = i(r + +r_) = ^stah 2 (^). (46) 

Using (|45p and (|46p. we can combine the surface conservation laws for individual ions (|36p and write out the 
electrochemical potential in terms of c* and <j>* to obtain 

e^f = eV s -{qV s \ogc* + z+wV s (j)*)~ z + c*^- (47) 
Ot on 

3w dc* 
e-^ = eV s • (wV s logc* + z + q\7 s <j)*) - — . (48) 



Since we have explicit expressions for q and w, (|47| and (|48[) are completely specified in terms of bulk field 
variables yielding surface conservation laws that can serve as effective boundary conditions for the bulk 
equations written in terms of these field variables. 



6.2 Concentrated Electrolytes 

Following Iglic and Kralj-Iglic (43j [44] and Borukhov et al. [45l |46l [47], we can write a free energy that 
accounts for steric interactions: 

F = U -TS 

5^ log <*+ ^1-^]^ log^l-^Cij +^z jCl 0-e 2 |V(/)| 2 ^ da. (49) 

As for concentrated neutral solutions, steric effects are included by adding a contribution from the solvent 
entropy (second term in integrand) and we have assumed that the molecules of all species present in the 
system (including solvent molecules) are of the same size. It is important to mention that using (f4"9"]l as the 
free energy is only one way of including steric effects. There are more sophisticated (and accurate) theories 
involving statistical density functional theory 48, 49, 50] [ST]. Unfortunately, these theories are generally more 
cumbersome to analyze and often require advanced numerical methods in order to gain physical insight 52J. 
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Next, we derive electrochemical potentials for the i-th ionic species following Kilic et al. [53 by taking 
the functional derivative of the free energy with respect to the z-th concentration: 



Mi = lo S 



i-E 



+ Zi 



(50) 



3 3 . 



This expression for the electrochemical potential leads to a modified Nernst-Planck equation for ion trans- 
port [S3] 



dc 



(51) 



Notice that, unlike dilute electrolytes, the transport of different ionic species are directly coupled via the last 
term in the flux. 

Following our previous procedure, we can derive the concentration profile of each ionic species within the 
double layer by using the fact that each of the fit must be constant in the normal direction. This observation 
leads to a linear system of equations for the ionic concentrations, which is easily solved to yield 



Ci{z) 



cxp m - Zi(p{z f 



(52) 



As for concentrated neutral solutions, the concentration profiles have a Fermi-Dirac-like form. 

A little further progress can be made for the special case of symmetric, binary electrolytes. The excess 
concentrations of the ionic species in this case are given by 



7± = c± 



=f z . lb 

e + 



1 + 2^ sinh 2 (z + i>/2 



(53) 



where c*, ip and v are as defined in our discussion of dilute electrolytes and concentrated neutral solutions. 
Unfortunately, integration of these excess concentrations for the concentrated electrolyte model does not 
yield a particularly simple expression for the surface excess concentration |52j : 



r± = T sgn(C)z + eW-log(l + 2^sinh 2 (z + C/2)) 



+ (l-v)e 



z + c coshu - 1 
1 + 2^ sinh 2 u 



1-1/2 



log (1 + 2^ sinh u) 



du. 



(54) 



Combining this expression for the surface excess concentration and the electrochemical potential for concen- 
trated electrolytes (|50p yields the surface conservation law for concentrated, symmetric binary electrolytes. 



6.3 Dominant Terms in Electrochemical Surface Conservation Laws 

In general, the relative importance of the terms in a set of surface conservation laws (written in any form) may 
depend on the choice of time scales and the magnitude of surface transport. In the context of electrochemical 
transport, we find that the boundary conditions applied in many theoretical studies are merely the leading 
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order form (in different asymptotic and physical regimes) of the surface conservation laws derived in the 
previous section. For instance, in induced charge electro-osmosis problems in dilute solutions at weak applied 
electric fields (TH1 HO], e<7 and ew remains O(e) quantities. Thus, compared to the normal flux term, the 
surface flux term is negligible and the time-dependent term is only important at short times, t — 0(e). In 
this situation, the double layer charging equation (|T7|) becomes [2H1 US] : 

Yt =a ^ (55) 

where a — c* is the bulk conductivity of the solution and time has been rescaled using t = t/e so that the 
dynamics are on the RC time scale [20] [19]. At t = 0(1) time scales, only the normal flux remains an 
0(1) quantity, so we are left with a "insulator" boundary condition for the electric potential: 

At higher applied fields or for highly charged particles [4] [6] [7] [5] , the dilute solution electrolyte model 
leads to surface excess concentrations as large as 0(l/e) for some of the ionic species in the double layer. As 
a result, surface currents becomes important and evolution of ionic concentrations within the double layer 
occurs on an O(l) time scale. In this situation, no terms in (|36p are negligible, so we must retain all terms 
in (J37J) and 



7 Summary 

In this article, we have presented a general formulation and derivation of surface conservation laws at 
microscopically diffuse interfaces. Because surface conservation laws form the crucial connection between 
bulk and interfacial dynamics, it is important to know that they are theoretically well-founded. Our work 
fills this apparent void in the literature and provides a solid theoretical foundation based on techniques 
from asymptotic analysis. Our analysis has also led to explicit formulae for the surface and normal fluxes 
involved in surface conservation laws. In addition to formulae for arbitrary interfacial models, we have 
presented a formulation of surface conservation laws for the important class of interfacial models derived 
using the principles of non-equilibrium thermodynamics. This formulation provides a method for developing 
interfacial models in a systematic and physically sound manner. Finally, we have demonstrated the derivation 
of surface conservation laws in the specific contexts of diffusive and electrochemical transport. We emphasize, 
however, that surface conservation laws are very general and apply to a wide-range transport processes. The 
basic approach is to first specify (or derive) a model for the interfacial layer and then to compute the surface 
excess concentrations and surface flux using the general formulae presented in this article. 
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Notes 



1 Notc that differential operators will also need their units restored. 

2 Cmax = 1/a 3 , where a 3 is the characteristic volume of each solute/solvent particle. 

3 This assumption may be relaxed in situations where there are processes intrinsic to the interface (e.g., chemical reactions). 

4 In dimensional form, mobility coefficients are related to diffusion coefficients via the Einstein relation: Dij = kTLij. With 
our choice of energy scale, the dimensionless Einstein relation is Dij = Lij. 

5 This expression for the chemical potential implicitly assumes that the system is locally in thermal equilibrium. Given 
this assumption, it is easy relate (125 I t to the definition from equilibrium thermodynamics by recognizing that the functional 
derivative of the total free energy with respect to concentration is merely the partial derivative of the free energy density with 
respect to concentration, which yields the local chemical potential of the system at each point in space. 

6 In dimensional form, Lij would be a scalar multiple of the identity matrix. 

7 The factor of 1/2 in the definition of q is present for mathematical convenience. The total surface charge density is 2eq. 
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